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Abstract 

A method to extract primary 7-ray spectra from particle-7 coincidences at exci- 
tation energies up to the neutron binding energy is described. Prom these spectra, 
the level density and 7-ray strength function can be determined. From the level 
density, several thermodynamical quantities are obtained within the microcanonical 
and canonical ensemble. Also models for the 7-ray strength function are discussed. 

1 Introduction 

The Oslo Cyclotron Group has established a method to deduce experimental level densities 
and 7-ray strength functions |1|. These data contain essential information on nuclear 
structure and thermal and electromagnetic properties. In the last couple of years several 
fruitful applications of the method have been reported [2-10]. 

The most efficient way to create entropy in atomic nuclei is to break J = nucleon 
Cooper pairs in the core. A beautiful manifestation of pairbreaking is the backbending 
phenomena in rapidly rotating nuclei. In this work, the pairbreaking process is studied as 
function of intrinsic excitation energy (or temperature). Of great interest is the quenching 
of the pair correlations, a topic which is directly connected to the nuclear level density. 

Important applications of nuclear level densities and 7-ray strength functions are the 
determination of nuclear reaction cross sections from Hauser-Feshbach type of calcula- 
tions. These cross sections are used as input parameters in large network calculations 
of stellar evolution, and in the simulation of accelerator-driven transmutation of nuclear 
waste. 



2 Experimental method and techniques 

The experiments were carried out using the ( 3 He,a7) and ( 3 He, 3 He'7) reactions at the 
Oslo Cyclotron Laboratory with a 45 MeV 3 He beam. Self-supporting metallic targets of 
rare earth nuclei were isotopically enriched to ~ 95% and had thicknesses of ~ 2 mg/cm 2 . 

The charged particles and 7 rays were recorded with the detector array CACTUS, 
which contains eight particle telescopes and 27 Nal 7-ray detectors. Each telescope is 
placed at an angle of 45° relative to the beam axis, and comprises one Si front and one 
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Figure 1: Comparison between observed primary 7 spectra for the 162 Dy( 3 He, 3 He'7) 162 Dy 
reaction (data points) and fits (lines) according to Eq. ([!]). 



Si (Li) back detector with thickness 140 and 3000 /zm, respectively. The Nal 7-detector 
array, having a resolution of ~ 6 % at E 1 = 1 MeV and a total efficiency of ~ 15 %, 
surrounds the target and particle detectors. In addition, two Ge detectors were used to 
monitor the spin distribution and the selectivity of the reactions. 

In order to determine the true 7-energy distribution, the 7 spectra are corrected for 



the response of the Nal detectors with the unfolding procedure of Ref. JTTJ. In addition 
random coincidences are subtracted from the 7 spectra. The set of unfolded 7 spectra are 
organized in a (E, E 7 ) matrix, where the initial excitation energies E are determined by 
means of reaction kinematics utilizing the energy of the ejectile. This matrix comprises 
the 7-energy distribution of the total 7 cascade. 

The primary 7 matrix can now be found according to the subtraction technique of 
Ref. JT2| , see data points of Fig. |TJ. The procedure is based on the assumption that 
the decay properties of the particular reaction-selected distribution of states within each 
energy bin are independent on whether the respective ensembles of states are directly 
populated through the nuclear reaction or by 7-decay from higher lying states. 

The idea is now to find two functions, the level density p(E) and the 7-energy depen- 
dent function F(E 7 ), that according to the Brink- Axel hypothesis [[L3], should describe 
the primary 7-ray spectrum at excitation energy E: 

P(E, E 7 ) oc p(E — E 7 )F(E 7 ). (1) 
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Figure 2: Observed 7-ray strength functions /ei+mi and level densities p. The data may 
be downloaded from |http: / /www. fys.uio.no/kjerne/english/ cyclo/compilation.html . 



Here, P{E 1 E 1 ) is fitted to the observed primary 7-ray matrix Q. In Fig. [l] the best fit 
to P is shown for the 162 Dy nucleus. All primary 7-ray spectra are seen to be very well 
described by the same p and F functions. However, from p and F one can construct other 
functions which give identical fits to the data by 

p(E-E 7 ) = Aexp[a(E - Ej]p(E - E 7 ), (2) 
F(E y ) = Bexp(aR r )F(E y ). (3) 

Here, the parameters A and a can be determined by fitting the level density to the 
number of known discrete levels at low excitation energy and to the level density estimated 
from neutron-resonance spacing data at high excitation energy. Furthermore, for dipole 
radiation, F is proportional to fE 7 3 , where / is the 7-ray strength function. Thus, 
according to Ref. |9[], the parameter B can be determined from the known total 7 widths 
of neutron resonances. Finally, the normalized level densities and 7-ray strength functions 
are shown in Fig. 0. 
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Figure 3: Left panel: The entropy as a function of excitation energy in 166 Er and 167 Er. 
Right panel: The temperature and heat capacity from the micro-canonical (data points) 
and canonical ensemble (lines). 

3 Thermodynamic properties 

High nuclear level density is due to the many ways the potential and kinetic energies 
of individual nucleons can add up to a certain excitation energy. Each valence nucleon, 
i.e. nucleons not coupled in Cooper pairs, carries a single-particle entropy of 1.5 — 2.0 in 
units of the Boltzmann constant ks- Thus, a broken Cooper pair contributes significantly 
to the total entropy and heat capacity. Observed structures in the caloric curve and heat 
capacity may therefore indicate the breaking of Cooper pairs. 

The microcanonical partition function is simply the multiplicity of nuclear states, 
which experimentally corresponds to the level density of accessible states. Thus, the 
experimental level density p(E) is our starting point for the extraction of thermodynamic 
properties of nuclei. 

The entropy in the micro-canonical ensemble is determined by 

S(E) = Hp/p ), (4) 

where the Boltzmann constant is set to unity (ks = 1), and po is adjusted to obtain S ~ 
in the ground-state band of even-even nuclei. The nuclear temperature and heat capacity 
is defined by 

T(E) = (f§) ' and (5) 

C V (E) = (§)". (6) 

Figure |3| shows the entropy, temperature and heat capacity for 166 Er and 167 Er deduced 
in the microcanonical ensemble. The small bumps in the entropy curves are seen to be 



4 



30 








^ 20 




\ ,/ 




J 10 


- 


/ ,7, Yb 


,,f ,7, Yb 


>> 

' 






/ 


8 40 
Cl 

8 30 








o 








CD nn 








10 




/' ,72 Yb 


,72 Yb 










" 


0.0 0.2 


0.4 0.6 0.8 


2 4 6 8 



T [MeVl <E> [MeV] 

Figure 4: Canonical heat capacity as function of T and (E(T)). 

enhanced through the differentiations performed in Eq. (]5|). The spectacular feature of 
negative branches in the heat capacity is a direct consequence of negative slopes in the 
caloric curve. Until today, their interpretation is a controversial topic. 
The canonical ensemble theory is based on the partition function 

oo 

Z(T) = £ u(E)e- E / T , (7) 

E=0 

which is a Laplace transform of the multiplicity of states uj{E) = AEp(E), where AE is 
the energy bin. The thermal average of the excitation energy in the canonical ensemble 
is 

oo 

(E(T)) = £ Eu{E)e- E ' T . (8) 

E=0 

By the Laplace transform in Eq. (|7|) much of the information contained in the micro- 
canonical level density is smeared out. Thus, fine structure in the thermodynamic ob- 
servables in the microcanonical ensemble will not be visible in the canonical ensemble. 
The lines in Fig. [5] display the smooth dependence of the canonical temperature on (E). 
The corresponding heat capacity is the derivative by 

Cv(T) = (9) 

and is shown in Fig. || for 171 > 172 Yb as a function of T and (E). The heat capacities show 
a pronounced peak as function of temperature at T c = 0.5 MeV, which we interpret as 
the critical temperature for the pair breaking transition. 

4 Gamma-ray strength functions 

In Fig. |5] the experimental 7-strength function is fitted by a theoretical strength function 
taking into account both the giant electric dipole resonance and the spin-flip resonance. 
In addition, a weaker resonance at lower energies is needed in order to fit the experimental 
data. Because of the much lower strength of this resonance compared to the giant electric 
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Figure 5: The experimental El + Ml 7-ray strength function (data points) of 166 Er (left) 
and 167 Er (right). The solid line is the fit to the data by the theoretical model. The 
dashed lines are the respective contributions of the GEDR, the GMDR, and the pygmy 
resonance to the total theoretical strength function. 



dipole resonance (GEDR), it is denoted pygmy resonance. The El radiation is described 
by UJ 



/El( ^ } " stt^v E^-E^y ' (10) 

where we apply T as a constant fit parameter. The Ml radiation is described by the 
Lorentzian 

where <tmi, Tmi and Emi are giant magnetic dipole resonance (GMDR) parameters ||. 
Furthermore, the pygmy resonance is described with a Lorentzian function f py (similar 
as Eq. (0)), where the pygmy- resonance strength o" py , width r py and centroid E py have 
been fitted in order to adjust the total theoretical strength function to the experimental 
data. The pygmy resonance parameters obtained for several rare earth nuclei are shown 
in Fig. |. The systematics indicate that the strengths a py of these soft dipole resonances 
are quenched when approaching the N = 82 shell gap. 

The pygmy resonance has been explained by the enhancement of the El 7-strength 
function |I6[ . Still, Ml character cannot be excluded. At an excitation energy around 
3 MeV, there is a concentration of orbital Ml strength in the weakly collective scissors 
mode |T7|]. This mode was first observed in electron-scattering experiments ||18|| , and is 



confirmed by the (7,7') reaction [19|. Clarification of the electromagnetic character of 
the soft dipole resonance observed in our work awaits new experimental results. 
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Figure 6: Systematics of the pygmy resonance parameters for odd (filled squares) and 
even (open squares) nuclei as a function of neutron number N. The resonance energy 
E py , the width T py and the cross sections <r py are shown in the upper, middle and lower 
panels, respectively. 



5 Conclusion 

Thermodynamic observables have been deduced from the level density and display sig- 
natures of phase-like transitions within the microcanonical and the canonical ensemble, 
interpreted as the transition from a strongly pair-correlated phase to an uncorrelated 
phase. In the microcanonical ensemble one may observe details of the successive breaking 
of nucleon pairs, information which is hidden in the canonical approach. The canonical 
ensemble on the other hand, reveals the average properties of the pairing transition. Us- 
ing the canonical Cy(T) quantity as thermometer, a local maximum is found at T c ~ 0.5 
MeV, indicating the breaking of Cooper pairs and quenching of pair correlations. 

The experimental 7-strength function is fitted by a theoretical strength function, as- 
suming that the 7 decay in the continuum is governed by dipole transitions. The contri- 
bution of electric and magnetic dipole radiation to the 7-strength function is recognized, 
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and a bump is observed in the 7-strength function at E 1 ~ 3 MeV. A measurement 
of the electromagnetic character of this soft dipole resonance will be crucial in order to 
understand the physics behind this exciting phenomenon. 
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